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ABSTRACT 


This study establishes the existence of damping structu- 
ral properties in the San Andreas Fault zone through an input 
model of surface loading. Any changes in the density, mass 
or bulk elastic properties of a fault zone will not only pro- 
duce a change in the gravity field of the earth, but also a 
change in the earth's response to the tidal field generated 
by the moon. Knowledge of these changes is important to en- 
gineering constructions, pre-programmed defense systems, etc. 
To adequately measure these responses, sites on opposite 
sides of the San Andreas Fault zone were selected and tidal 
gravity was continuously measured for finite periods of time. 
The Pacific Ocean and its very low frequency lunar portion of 
the ocean load spectrum was used as a driving function to 
compare observed data with theoretically generated data for 
the same sites. A least squares criteria was used to fit the 
observed and theoretical data to the nine most important 
spectral lines and to obtain the phase and amplitude for each 
line at both sites. Utilizing the most stable lines to com- 
pare values, a significant drop (about 10%) in both phase 
and amplitude occurred between the two sites. The anomalous 
decrease in both areas of analysis discloses a structural 


damping in the San Andreas Fault zone. 
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i. INTRODUCTION 


A. STATEMENT OF THE PROBLEM 

The combined gravitational fields of the moon and the 
sun acting on the solid earth are sufficient to displace its 
surface on the order of 20 cm. semi-~daily fi a Untortunate— 
ly, this direct forcing function is so broad in its areal 
event (1.e., Several thousand kilometers) that it cannot 
reasonably be used to study responses of small crustal 
structures (such as fault zones) [2] which are on the order 
of a few kilometers in width and depth. However, the response 
of the ocean to the lunar~-solar gravitational field (on the 
order of meters) constitutes a sufficiently energetic forc- 
ing function such that under certain conditions dynamic 
properties of crustal structures can be studied. 

The conditions under which the ocean can effectively be 
used as a forcing function to study the properties of a 
Shallow structure are defined as: 

1. The distance between the edge of the ocean load and 
the fault structure should not exceed the smallest dimension 
of the fault zone by one order of magnitude [3]; 

2. Noise response in the earth tide records obtained 
during this experiment must be kept to a minimum. Minimiz- 
ing noise response is achieved by using parts of the ocean 
tidal spectrum which are both stationary and relatively un- 


affected by meteorological influences. 











The two lines chosen for amplitude and phase analysis in 
this study are the principal lunar diurnal 07 and the princi- 
pal lunar semi-diurnal M, lines which are the most energetic 
in the theoretical tidal spectrum and are the most removed 
from possible contamination due to meteorologically induced 
MewsesS. 

Earthquake prediction research explores the causes of 
earthquakes from two points of view: (1) failure brought 
about by an external increase of shear stress which exceeds 
mee shear modulus of the fault; (2) failure brought about 
when the internal properties of the fault change (coefficient 
of shear) until they are reduced below that of the external 
shear being applied. 

in this study responses to tidal gravity were measured 
at two points on either side of the San Andreas Fault. Tidal 
pravity is a vertical measure of the acceleration due to 
gravity. The value of gravity is on the order of 10? em/sec/~ 
sec or 10° gal. Resolution on the order of 107° gal is 
achieved with the instruments used (described below) ee 
mesultsS accurate to one part in a 107, The distance between 
the two stations (about 5 kilometers) is small compared to 
the distance to the edge of the ocean load (Pacific Ocean, 
about 30 km). The amplitude and, phase of the two measured 
tena tidal constituents (0, and M, ) were compared with the 
amplitude and phase of the same two lines derived from an 


analysis of the theoretical tidal gravity calculated at each 








of the two sites for identical periods of time. The differ-— 
ence in the observed phases between the sites was compared 
With the difference in theoretical phases between the sites 
and a 10% change in the means was observed. 

This result suggests that the principal source for the 
phase difference between the observed and theoretical values 
at both sites is due to the response of these sites to ocean 
tidal loading at the same frequency. However, a 10% devia- 
tion of the difference between these two sites from the mean 
is significantly large and is not due to a simple shift ofa 
few kilometers further away from the coastline. Instead, the 
anomalously large difference in the comparison is explained 
by the fact that the stations are not just separated by a 
few kilometers but that they lie on opposite sides of a weak-— 
ness in the earth's crust. The absolute difference between 
the most inland station and the other is a lag of about 
three degrees at the M, frequency. This is interprted pri- 
marily as a viscoelastic response of the fault zone to ocean 
tidal loading. The inland edge of the San Andreas Fault 
reaches it peak amplitude of ocean tidal response approxi- 
mately 1/10 hour after the peak response of the ocean side 
of the fault. 

While this study did not attempt to directly calculate 
any geological strength parameters, it did suggest the 
existence of a crustal weakness with viscoelastic properties 


that strongly dampen the M, line er the occan tidal boading. 
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B. BACKGROUND 

There exists in any multi-body system a measurable gravi- 
tational potential due to the attacting fields of the many 
masses (or bodies) present. Since all bodies have motion 
relative to all others ina system, the potential fields 
fluctuate with time. On the earth, only the sun and the moon 
contain sufficiently large quantities of mass at a close 
enough distance to significantly affect tidal variations or 
time aie eteions nieuwe PpOuenulal field. 

The effects from the sun and the moon generate earth 
tides in a similar manner to their production of ocean tides. 
Unfortunately, the study of earth tides is complicated by 
the presence of oceans which also cause changes in earth 
tides. Although quite small, ocean loading effects on the 
earth tides have been detected in measurements as far inland 
as mid-continent in the United States [4]. The mass of the 
oceans complicates measurements of the earth tices in the 
same ways as external masses (the sun and moon) plus (1) 
the loading effect due to water (additional mass) and (2) 
deformation of the earth due to the redistribution of ocean 
waters as’the water more freely fluctuates. 

The combination of effects of ane sun, moon and oceans 
Over time results in a continuous change in the earth's 
gravity. The gravitational force at any given point in time 
and space can be broken into two components, one normal to 


the plane tangent to the earth's surface and one in the 


at 





plane tangent to the surface. This study will concentrate 
on the former component called tidal gravity. 

Past geological and seismological studies have examined 
the bulk elastic properties of the whole earth and the dis- 
tribution of these properties in the core and the mantle. 
Few studies have narrowed their objective to examining the 
bulk elastic properties of only the crust of the earth and 
in a relatively small area. Seismic studies show that the 
San Andreas Fault zone has brittle properties to a depth of 
only ten kilometers, which is neither the bottom of the 
falc nor the bottom of the continental plate [5]. It is 
approximately one-third of the depth of the plate. Authors 
of these studies quite often generalize their results to 
large area averages and tend to overshadow the relatively 
smell area of a shallow fault zone in the final analysis [2]. 
The seismic studies and their related high frequencies (from 
earthquakes) also do not allow determination of anelastic 
responses, but very low frequency studies such as those re- 
sulting from ocean loading should disclose any anelastic 
properties present. 

This study is a portion of a larger project, (headed by 
Dr. M. D. Wood) at the National Center for Farthquake Re- 
search (NCER), to establish that the response of the earth 
to tidal loading is significantly perturbed by the San 
Andreas Fault and also to determine through spatial and 


temporal means the bulk properties of the fault zone. 
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There are two possible sources of error in analyzing 
these data that should be noted. First, due to many con- 
straints, the three experiments in this study were performed 
sequentially instead of simultaneously. Second, the San 
Francisco Bay, due to its size and its proximity to the sites, 
will affect the values measured in the experiment. The ef-— 
fects of these errors are assumed to be relatively small, 


but their exact effects are undetermined. | 
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TI. DATA ACQUISITION AND PROCESSING 


A. LOCATION SELECTION 

The locations selected for the two main sites for this 
analysis were obtained from a threefold criterion. First, 
it was desired to ensure that the two sites were on opposite 
sides of the San Andreas Fault zone and nearly equidistans 
from the fault region. This was necessary in order to ob- 
tain spatial analysis of the tidal wave as it passed from one 
Site to the other by traversing through the fault region. 
mwecond, the sites were desired to be on a line nearly per- 
pendicular to the coastline. Third was the availability of 
suitable sites. It was imperative that the physical location 
of each instrument be based on a very sturdy foundation, 
ideally on deep lying rock, to minimize noise. To further 
eliminate noises (in the form of other energies such as wind, 
temperature fluctuation, surface vibrations, etc.) it was 
dost. age to place each recording instrument below the sur- 
face of the earth. Time and financial restrictions forced a 
compromise of this desire to a level, surface (or very 
shallow) base that was protected from as many of the other 
noise—producing sources as possible. In addition to these 
fear ictuions, the security of the expensive recording instru- 
ments became a major concern. The instruments could not be 
left in the open where vandals or sightseers might tamper 
with them. It was also desired to have daily time marks 


entered on each record, generating the requirement for 


1h 











Ew@ieene FG be present at the site each day to perform this 
task. Appropriate sites with the necessary attendance re- 
quirements were located [Figure 1]. One recording instru- 
ment was in a basement vault on a base on concrete poured on 
solid rock (Block experiment); another was in a separate 
garage on a concrete base poured on rock (Rose experiment). 
Both properties were private and actual instrument locations 
Were secure. The third site was located in the NCER Build- 


ing in Menlo Park on a concrete floor (NCER experiment). 


B. EARTH TIDE RECORDINGS 

All the earth tide recordings were taken from LaCoste— 
Romberg Earth Tide Gravity Meter type instruments [1] on 
paper strip charts [Figure 2]. These were recorded from 
NCER owned instruments run by 120v AC with 12v automobile 
batteries as a back-up. 

wovallow iuture adjustments for inaccuracies in chart 
Eo sad , stretching or shrinking, etc., daily time marks were 
mechanically entered on each chart and labeled. Chart speed 
for each experiment was two inches per hour with a displace- 
ment of one inch equal to 30 yrgal. Each experiment was 
attempted to run at least thirty days to give a large enough 


sample for adequate analysis [6]. 


CGC. DATA PROCESSING 
The strip charts were digitized onto magnetic tape by 


mecwagneamcalma 502 digitizer. The cross-—hair cursor moved 
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vertically or horizontally as desired across the strip 

chart taking discrete samples every .Ol inches and recording 
the readings. A special indication was noted at each daily 
time mark entered on the chart by monitoring personnel, and 
the time and date were recorded by the person performing the 
digitization. The digital tape data was later corrected 
momechar. Stretching, Shrinkage, etce., by proportionally 
distributing the actual recorded values between time marks 
to the correct number of .0O1 inch increments (.01 inch 
corresponds to 18 seconds of time) required between the daily 
time marks. This localized physical errors between time 
marks (usually every 24 hours) and minimized accumulation 
@ieerror. 

The digitized data consisted of readings based on a 
continuous time record. These data were then processed with 
a computer program (designed by M. D. Wood and P. A. Marshall 
[7] from the National Center for Earthquake Research (NCER) 
and modified by M. D. Wood and R. L. Keller for use in this 
study) which decimated the sample interval from .01 inches 
to hourly values taken on the hour. 

sampling theory requires a sample interval at least twice 
as often as the highest expected frequency. Implicit in this 
is the assumption that the recorded data be of a reasonable 
duration of time. A reasonable time duration is propor- 
tional to the reciprocal of the smallest difference between 


any two frequencies expected in the data [6]. Since only 
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Nine lines of the 80-plus theoretical lines are used, it can 
be calculated that at least one month of data is necessary 
for a 1% resolution of tidal data [6]. The smallest fre- 
quency line used in this study is the M), with a frequency 
of less than 4 cycles/day, requiring at least & samples per 


day. By oversampling, the resolution remains well within 1%. 


IIT. LEAST SQUARES ANALYSIS 


A. AMPLITUDE ANALYSIS 

The observed tidal gravity records were fit with a least 
squares computer program to the nine theoretical lines of the 
earth tide spectrum (explained in Appendix A). The computer 
program (designed by M. D. Wood [6] and modified by M. D. 
Wood and N. E. King of NCER and R. L. Keller) performs a re- 
gression of the data to fit both an neh polynomial anda 
Sine wave function simultaneously. The results give the 
amplitudes and phases of the nine lines mentioned above. 

It is a well established fact that the principal lunar 
semi-diurnal M. and the principal lunar diurnal 05 lines of 
the earth tide frequency spectrum are the most invariant of 
feecne tidal lines [1], [6], [8]. Of these two, the M, 
jine has proven to be less variant from its theoretical value 
than the other lines in previous studies conducted in this 
Perion (6]|, [7], [8] and was thus used as a normalizing fac~ 
tor for all of the spectral lines. The normalized amplitudes 


were computed using the following formula: 


aR 








| & 


A; = As/Aye 


where Ay is the (dimensionless) relative amplitude, A. is 
the amplitude of the constituent line being normalized, and 
Avo is the amplitude of the M, line. This normalizing pro- 
Céss results in dimensionless values, facilitating compari- 
son between observed and theoretical results. 

The normalizing of the amplitudes of the selected 
lines in the earth tide spectrum was conducted for both the 
imeoretical data |Table III] and for the observed data [Table 
IV]. Visual displays of these data are given in Figure 7 
and Figure &. 

Taking the expected amplitude (theoretical) of each 
constituent line for each of the sites and subtracting the 
Sormesponding observed amplitude, the amplitude differences 
are obtained [Table I]. Differences can be expected due to 
forcing nine lines of a frequency spectrum to absorb what is 
acknowledged to consist of at least 89 lines and also due to 
other non-tidal factors such as seasonal noise, extreme wea- 
TMer conditions, etc. These differences should remain 
about constant between sites (resulting in low values in the 
last column) and would indicate a high stability of a con- 
Stituent line. The last column of Table I shows a standard 


deviation of about 5% for the 0, line from its average theo- 


al 
retical mean. This indicates that the 0, line could be pre- 
mtewed to have a measured (observed) value with a standard 

deverwion equal to about 5% of the theoretical value of the 


amplitude. 
20 














The normalizing of the amplitudes of the selected 
lines in the earth tide spectrum was conducted for both the 
theoretical data [Table II] and for the observed data [Table 
III |. Visual displays of these data are given in Figure 3 
and Figure 4. These tables and figures are presented here 
to facilitate comparison of the results obtained in this 
study to those of other completed studies [6], [7], [8] 
making use of similar types of displays. 

A relative gravimetric factor te Can. be “CoMmpuped joy 


the following formula: 


Ce 


> AG (experimental)/A; (theoretical) 


where the Ay are described above. Using the two most stable 


lines to establish a comparative base, we compute: 


G" (0, /M, ) = A* (0, /My ) (experimental )/A* (0, /M, ) (theoretical) 


lm@earetically, this ratio should equal to 1.00, but with 
other dynamic considerations a 1% variation [9] is accepted. 
Ten world-wide stations studied by Harrison and others [10] 


resulted in an average value of : 
te 
G(0,/M,) = 0.989+0.035. 


Four stations in the San Francisco Bay area analyzed by Wood 


[6] gave an average value of: 


G"(0,/M,) = 0.938+0.026. 


Zi 











The average value obtained by the three stations in this 


study gave: 


G"(0,/M,) = 1.05110.061. 


These values should not be different, regardless of loca- 
tion, yet this study and that by Wood indicate values about 
h% to 5% outside of the accepted range. Both of these 
studies were conducted on the West Coast near the San 
Francisco Bay Region and could indicate a local disparity 
from areas with thicker crustal depths (The California coast- 
al regions are known to have a much thinner crustal depth 
than most other continental areas). However, no conclusion 
is drawn from this study other than to mention that this dis- 


parity should not exist. 


B. PHASE ANALYSIS 

The observed tidal gravity records were simultaneously 
fit in the least squares program (described in the previous 
section) to the nine theoretical waves of the earth tide 
Spectrum to estimate the phase of each constituent line. 
The values for both theoretical and observed phases for the 
M, Jines are given in Table IV. All charts were maintained 
in reference to Greenwich Time (GMT) except for the Rose 
Site which was marked in local daylight time. A seven hour 
difference existed between these base reference times and 


the phase was appropriately advanced for this discrepancy. 


“ie 





The results show a phase difference between the theoreti- 
cal and observed Me line at both the Roge site and the Block 


site. 
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Figure 3. Composite Theoretical Amplitude Histogram 
Normalized by M, Component 
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Figure 4. Composite Observed Amplitude Histogram, 
Normalized by M, Component. 
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The average value of G" (0, /M, ) was much higher than those 
produced by the experiments done by Harrison and others [10] 
and Wood [7]. As noted above, no explanation for this has 
to date been found. The relative gravimetric factor showed 
a downward change between the two sites. It decreased by 
more than 0.1, or about a 10% drop in value. 

In the examination of the two most Ste Ne lines in the 
earth tide spectrum, the value of the 0, /My (the normalized 
amplitude of the O05 constituent line) decreased from the Rose 
to the Block station by a value of .0488 which is about a 
7e% decrease between the two sites. Theoretically, the nor- 
malized 0, line amplitude should not show much change in 
value, specifically across a distance of 5 km. From the 
values in Table II, it can be shown that about a 2% increase 
should have occurred. 

There also exists a phase difference between the theore- 
tical and observed data at each site. The observed and 
theoretical lines should all be in phase with each other. 
Allowing for small variations due to noise and tectonic in- 
terferences prior to the arrival of the M, line at the Rose 
site (again the most stable line is selected for comparison), 
there may exist a phase difference between theoretical and 
observed lines. Between the Rose and Block sites, however, 
this difference should remain constant. In fact, it does 


MOGs) 1 changes by more than 23 degrees. This change is a 


owt 





decrease in phase difference, or a lag, of about 93% between 
the two sites. Physically interpreted, this is a time delay 
of about 1/10 hour for peak amplitude of the M, wave to 


travel from the Rose to the Block site. 





BZ 


V. CONCLUSIONS AND RECOMMENDATIONS 


With the evidence advanced in this study, it appears 
that some disturbing factor has caused a large percentage 
decrease in tidal amplitude and phase between the Rose and 
the Block sites. These percentage changes are much too great 
to be explained by the 5 km. distance between the sites. The 
fact that the San Andreas Fault is located between these two 
sites leads one to suspect that drastic changes in earth 
tides between the two sites might be due to the fauit region 
and its different structural properties from surrounding area 
Properties. 

Assuming that the fault region is responsible for these 


tidal changes, it is possible to calculate the approximate 





strength parameters of the geological structure of the fault 
region itself through a more detailed geological study of 
Sites across the fault. 

since this experiment was limited in time, money, equip- 
ment, available sites, etc., it cannot provide absolute 
proof that the fault region is responsibie for the observed 
changes. The study provides evidence that the fault might be 
reason for the observed changes. It is recommended that 
further investigations of a similar nature be conducted at 
pmees VOUQ Of an inland bay to verify the results of this 
study. These future studies could also be set up to simul- 
taneously examine the specific geological properties exhibited 


by the fault region in each area of study. 
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APPENDIX A 
THEORY OF EARTH TIDES 


A study of the earth tide curve shown in Figure 5 reveals 
the presence of at least two frequency components of the 
earth tide spectrum. Identified in this figure are the diur- 
nal frequency with a period of approximately 28 days and the 
semi-diurnal frequency with an approximate 14 day period. A 
Fourier analysis of a record such as this reveals that there 
exist many constituent waves that combine to give us the wave 
form. Analysis of longer records would reveal semi-annual, 
annual, and even longer period waves which correspond to mo- 
tions of the sun and moon in relation to the earth. 

When the complex relationships between the moon (and sun) 
of a rigid earth are considered, the constituent waves (lines 
in the spectrum) derived via the Fourier analysis become 
Meemineful., Let us examine a fictitious model [6], [8] in 
which the moon is orbiting about the earth in the plane of 
the earth's equator at a constant distance with a constant 
velocity equal to the mean distance and velocity of the actu- 
al moon. The earth is also revolving on its axis once every 
24 hours and this motion plus the moon's gravitational at- 
traction on the earth's surface produces a tide with a per- 
iod equal to a lunar half-day and is described by the M, 
principal lunar semi-diurnal line, the most energetic of the 


semi-diurnal lines. In similar fashion, the So jonealalelaly estab 


3h 








solar semi-diurnal line is produced by a fictitious sun ina 
eircular orbit at a constant velocity and distance. 

If the moon's declination in its orbit (above or below 
the plane of the earth's equator) is now added to our simple 
model, the distance between the moon and an observation point 
on the earth varies with time, but ina different manner than 
the simpler model. In the simple model, the distance varia- 
tion with time was due only to the observer's position and 
the lunar position in the equatorial plane. Now a declina- 
falen Gistance away from the equatorial plane also varies with 
time as well as the consideration for the earlier model, add- 
ing another variable to be considered in the distance varia- 
tion between the moon and an observation point on earth. 

For observation points on the equator, the effects of the 
lunar orbit in the northern and southern hemisphere are equal 
in magnitude. For non-equatorial observers, the amplitude 
of the two diurnal maxima would vary in time with the chang-~ 
ing distance from the observation point to the moon. This 

is referred to as the diurnal inequality. The cyclical mo- 
tion of the moon's declination causes a fortnightly modula- 
tion of the earth tides. 

macame tne eclliptuacaty of the moon's orbit to our model 
now leads to the cyclical changes in the distance from tne 
moon to the earth and causes modulation of both the amplitudes 


and phases of all the above effects. 
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Associating a theoretical wave with each of the motions 
of the moon does not complete an explanation of the earth 
tide spectrum. Secondary waves are apparent due to the beat 
frequencies of the principal tide waves. These beat fre- 
quencies are obtained from the sum and difference frequencies 
of the two principal waves. When the difference in the lunar 
and solar principal semi-diurnal waves (M, and S,) is taken, 
a fortnightly modulation occurs at the beat frequency. The 
junar elliptic semi-diurnal wave, No, and the M, wave com— 
bine for a beat frequency to form a modulation with a monthly 
period that is associated with the eccentricity of the lunar 
Cieit. The Ko luni-solar semi-diurnal declinational wave is 
composed of lunar and solar components which are analytically 
inseparable, as is the Ky luni-solar diurnal declination wave. 

Further beat frequencies exist which yield waves derived 
in a similar fashion. Table I lists the principal theoreti- 
cal: tidal waves of the earth tide spectrum which are of 
importance to this study. The amplitude of each of these 


waves is dependent on latitude and varies with time. 
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APPENDIX B 
THEORETICAL TIDE GENERATION 


The gravitational potential at some point P on the sur- 
face of the earth [Figure 6] due to another mass (moon or 


pom) can be expressed as [6], [11}, [12]: 
CoO oO 
ez SE 3 
Ne =. 
W = ey at e ee Zz ) = Wi 
n=o n=o 


where G is a gravitational constant 

M 1s the external mass 

16 is the radius of the earth 

R is the absolute distance from the earth 
to the external mass 

ie 


(z) is the Legendre polynomial 


n 
Z is the angular position of the external 

mass relative to the zenith of point P 
is is a point on the earth's surface. 


Expanding the terms in the summation above necessitates the 
calculation of the Legendre polynomials. The first few terms 
of the Legendre polynomials are [12]: 

Piz) =1 

P, (2) 

Po(z) = $(32°-1) 

P3(z) = $(52°-3 ) 

2) a 1/8(352*-302°43) 


Z 


Examining these results it is seen that the first term, 
W =GM(1/R) , is constant and has no effect on changing the gravi- 
tational potential. The second term, W, =GM(r/R* )cos(z), uni- 
foumiy anteeus all points on the surface of the earth and has 


Monetiect on changing the potential. The third term (and all 


De 





remaining terms) affect the tidal potential and are ail 
candidates for significant effects in the summation. The 
r/R term in the summation is approximately 1/60 for the moon 
and 1/2400 for the sun and after n=3 accounts for less than 
2% of the previous term. Thus, the gravitational potential 


can be reduced to the simple expression: 
2 fe 
be r i Z) 
Wy = GM,3 Po(2) = GM,3 2(3cos (2)-1)- 
To express z in terms of fundamental angular parameters 


of position, spherical trigonometry is applied where: 


& = declination of external mass M 

2 howimean edhe somr pO lain 

® = zenith angle of external mass at point P 
ro lacitudesoul! 

A. = longitude of P 


Z= polar angle of external mass M 
VSM Caestrigonometricerelation: 

cos(«)= cos(e)cos(z) + sin(e)sin(z)cos(t) 
cos (YP) 


€o5( Zecca 2) = .cos (4 ) 


and cos(e)= sin(Y ), sin(e) 


the gravitational potential becomes 


Zc 
Wo=GM— 2 (2ieant) sin( 6 J costs )cos( Veo ny 
R 


end simplaties to 


Wy = 3K(sin®(?)-1/3)(sin*(6 )-1/3) 
PC sine sante > \eoslt )) 
- K(cos*(£)cos“(6 )cos(2t) ) 
where K = 3/4CMs 
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This equation is shown in Figure 7. 

Using the moon as an example, the first term in this ex- 
pression contains the variable & (fis constant for a given 
point P), namely Hines) TweLrcwetsoeani even Lunction varying 
as the declination of the mass M traverses from its maximum 
positive value to its maximum negative value and back each 
month. The functional value stays positive (since Sls) 
is an even function) and the frequency of- the first term be- 
comes two cycles per month and is called the fortnightly wave. 
The second term varies with § and t, but & completes a cycle 
in one month while t has a period of about one day. There- 
fore, this term has a frequency of one cycle per lunar day 
and is called the diurnal wave. Similarly, the third tern, 
with a frequency of two cycles per lunar day, is called the 
semi-diurnal wave. 

A similar situation occurs with the sun as the external 
mass, and the sum of the results from the moon and the sun 
yields the terms of the tidal potential. The remainder of 
the tidal potential spectrum consists of the higher ordered 
terms and combined lunar-solar terms neglected in this analy- 
Sis. For the purposes of this study the nine waves described 
in Table V are sufficient since they contain more than 90% 
of the effects of the 89 theoretical lines. 

In the generation of the theoretical earth tides for 
this study, J. C. Harrison's program [13] was used. In the 


development of his model for generating earth tides, Harrison 


bY 





uses Wy for the moon (as described above) and Wy for the sun 
and also includes We tor Chemmeon FCO achieve greater preci-— 


egom in his results. 
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Geometric Relation between Station Coordinates and 
Tide Producing Body Angular Coordinates. Wherex 
is the zenith angle of the tide producing body at 
the station location P. P is at colatitude™ and 
longitude A. Z is the polar angle and ? the 
terrestrial east longitude of M. (Adopted from 
Munk and Cartwright ) 
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